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Quantum Field Theory with fields as Operator Valued Distributions with adequate test functions, 
-the basis of Epstein-Glaser approach known now as Causal Perturbation Theory-, is recalled. Its 
recent revival is due to new developments in understanding its renormalization structure, which 
was a major and somehow fatal disease to its widespread use in the seventies. In keeping with the 
usual way of definition of integrals of differential forms, fields are defined through integrals over 
the whole manifold, which are given an atlas-independent meaning with the help of the partition 
of unity. Using such partition of unity test functions turns out to be the key to the fulfilment of 
the Poincare commutator algebra as well as to provide a direct Lorentz invariant scheme to the 
Epstein-Glaser extension procedure of singular distributions. These test functions also simplify the 
analysis of QFT behaviour both in the UV and IR domains, leaving only a finite renormalization 
at a point related to the arbitrary scale present in the test functions. Some well known UV and IR 
cases are examplified. Finally the possible implementation of Epstein-Glaser approach in light-front 
field theory is discussed, focussing on the intrinsic non-pertubative character of the initial light-cone 
interaction Hamiltonian and on the expected benefits of a divergence- free procedure with only finite 
RG-analysis on physical observables in the end. 

PACS numbers: ll.10.Ef, 11. 10. St, 11.30.Rd 



I. INTRODUCTION 

Causal perturbation theory (CPT) goes back to ideas 
from Stueckelberg, Rivier and Green [1], Bogoliubov, 
Parasiuk and Sirkov [2] and has been developed by 
Epstein and Glaser [3] and more recently by Scharf [4]. 
A little after Epstein and Glaser came formulations [5] 
to construct finite gauge theories, based on regulator 
free momentum substraction schemes. In the last 
publication the questions of the finiteness of Lorentz- 
invariant amplitudes and of symmetry conservation 
are treated seperately. Moreover, as opposed to CPT, 
the importance of causality issues in relation to the 
finiteness problem itself is not of main concern. In 
CPT the S-matrix is constructed as a formal functional 
expansion in an unimodular interaction-switching test 
function. Whereas in the traditional approach the 
S-matrix amplitudes are determined by equations of 
motion, in CPT they are determined inductively by 
causality conditions imposed on switching test functions. 
Such causality conditions turn out to be as stringent 
for the S'-matrix evolution as the equations of motion 
themselves. On an heuristic level it may be seen as 
follows: if gi(t) and gi (t) are two such switching test 
functions with supports such that supp[gi] £ (—00, s) 
and supp[g2\ € (s, +00) respectively, then causality 
implies that S(gi + g^) = S(gi) ■ S{g2), which is the 
functional equation for an exponential. 



The S'-matrix amplitudes are time ordered products of 
operator valued distributions (OPVD). They are split 
into causal advanced and retarded pieces, in such a way 
that any singular behaviour at equal space-time points 
is avoided. The procedure has long been recognized 
as mathematically rigorous [6, 7] since it is divergence 
free at each iteration stage. Nevertheless Epstein-Glaser 
work almost fell into oblivion among QFT practioners 
mainly because 

(1) in its original version there were difficulties in disen- 
tangling the multiplicative structure of renormalization 
and 

(2) the rise of renormalization group methods turned 
research interests into this direction. 



However Epstein-Glaser approach still enjoys popular- 
ity in the world of mathematical physicists dedicated to 
the construction of a rigorous and mathematically well 
defined QFT [8-10]. In this context the key issue is the 
extension of singular distributions to the whole space- 
time manifold. In mathematics a rigorous way to define 
an extension of a singular distribution is a weighted Tay- 
lor series surgery: one throws away an appropriate jet of 
the test function at the singularity and defines the ex- 
tended distribution by transposition at the level of the 
functional from the corresponding Taylor remainder of 
the test function. Transposed to Fourier space the pro- 
cedure amounts to a substraction method [10] which in- 
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eludes BPHZ rcnormalization [2, 11] as a special case. In 
a Minkowskian metric this is equivalent to the implemen- 
tation of causality while in its Euclidean counterpart it 
is a symmetry preserving prescription for substractions. 
In light-front quantum field theory (LFQFT) there is a 
compelling reason to introduce test functions related to 
the consistency of the canonical quantization scheme it- 
self. It is best seen for the massive scalar field. The 
form of the LF-Laplace operator leads to a hyperbolic 
equation of motion which requires initial data on two 
characteristics. Canonical quantization in terms of ini- 
tial field values in the light cone time is however possible 

provided [12] lim p +^ X p P + - = 0, where x(p + ) i s the field 
amplitude at p+ = p° +p 3 . With fields as OPVD this 
relation becomes lim p +^ - , which is surely satisfied 
for the class of test functions f(p + ) used for the Fock 
expansion of the field operators [13, 14]. 

Our earlier approach with test functions [13, 16] in fact 
implements Epstein-Glascr treatment of singular distri- 
butions in the LFQFT context. Here we want to as- 
sess our particular choice of test functions as 'partition 
of unity' and show the resulting simplifications brought 
about for the treatment of Poincare's invariance and for 
the Lorentz invariant Taylor series surgeries discussed in 
[10]. 

In Section II the general definition of fields as OPVD is 
given. Due to specific properties of Euclidean manifolds 
it is shown that the class of admissible test functions can 
be restricted to a 'partition of unity', in close analogy 
with the usual atlas- independent definition of integrals of 
differential forms. Thus a QFT construction is possible, 
which is test-function independent. As a consequence, 
for scalar field theory it is shown that Poincare's com- 
mutator algebra is fulfilled, in contrast with the case of 
arbitrary test-functions. In Section III the Bogoliubov- 
Epstein-Glaser expansion of the S'-matrix is recalled. It is 
given in terms of the interaction-switching test function 
g G V(M A ). The proper Epstein-Glaser time ordered 
construction of the expansion coefficients is described. 
The implementation of Lorentz covariance of the proce- 
dure is presented and the simplifications brought about 
with test functions vanishing at the origin with all their 
derivatives are shown. They concern Lorentz covariance, 
the Taylor surgery itself and its interpretation in terms 
of the BPHZ rcnormalisation scheme. In Section IV, the 
use of Lagrange's formula [8, 10] for Taylor's remainder in 
connection with partition of unity test functions is shown 
to simplify also the analysis of singular distributions both 
in the IR and UV domains. A direct application of the 
method in the IR is shown to be successful for the free 
massive scalar field theory developed in the (convention- 
ally hopeless) mass perturbation expansion. It should 
also prove relevant to many problems of massless confor- 
mal field theories (CFT). Finally, following [17], in the 
UV the alternative interpretation of the effects of parti- 
tion of unity test functions in terms of Pauli-Villars type 
subtractions at the level of propagators is pointed out 



in Euclidean and Minskowskian metrics. In Section V 
Epstein-Glaser separation method of Section III for dis- 
tributions into advanced and retarded pieces (splitting 
procedure) is rephrased for the special case of partition 
of unity test functions. Their use simplifies greatly the 
derivation of the relation between the splitting proce- 
dure and dispersion relations. In Section VI a general 
discussion is given on the possibility of applying BSEG's 
method in a non-perturbative Light Front Quantization 
framework and on the benefits one might expect from 
the procedure. Finally some conclusions are drawn in 
Section VII. 



II. FIELDS AS OPVD, PARTITION OF UNITY 
AND POINCARE COMMUTATOR ALGEBRA 

A. Fields as OPVD 

The history of fields as OPVD is almost as old as Quan- 
tum Field Theory itself. In fact the mathematical devel- 
opments of Distribution Theory (or Generalized Func- 
tions) [18, 20] in the 1960's were intimately related to 
problems raised at that time by pioneers in the QFT for- 
mulation. This mutual feedback has lead to QFT formu- 
lations [21] in the 1960 — 70's which aimed at dealing in a 
mathematically consistent way with intrinsic divergences 
crippling the usual QFT perturbative analysis. Among 
these approaches the work of Epstein and Glaser [3] has 
long been recognized as mathematically rigourous and 
free of undefined quantities. It is still generating an im- 
portant literature in mathematical physics and we shall 
summerise here some of its recent developments. 
To introduce fields as OPVD one may consider, with- 
out loss of generality, the free massive scalar field in D- 
dimensions. The general solution of the Klein-Gordon 
equation is a distribution, ie an OPVD, which defines a 
functional with respect to a test function p{x) belonging 
to V{R D ), the space of Schwartz test functions [18], 

*(p) =< (t>,p >= J d^y^y)p(y). (HI) 

Here &(p) is an operator- valued functional with the pos- 
sible interpretation of a more general functional <fr(x,p) 
evaluated at x = 0. Indeed the translated functional is a 
well defined object [18] such that 

T x $(p) = <T x (j),p>=<<f>,T_ x p> {11.2) 
= J d^y^y)p(x-y). (II.3) 

Now the test function p(x — y) has a well defined Fourier 
decomposition 

P{X y) = / ^ iq(x - y) fill I 2 )- (H.4) 
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It follows that 
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e~^S(p 2 -m 2 ) X (p)f(plp 2 ). 



Due to the properties of p, T x Q(p) obeys the KG-equation 
and is taken as the physical field <fi(x). After integration 
over po and with u> 2 
taken as the (D — 1)— Euclidean integral 



p 2 + m 2 the quantized form is 



<p(x) 



(2tt)(' D - 1 ) 2uj p 



a+e rpx +a p e' rpx }. {11.5) 



f(uj 2 ,p 2 ) acts as a regulator [13, 14] with very specific 
properties [15]. This expression for ip(x) is particularly 
useful to define correlation functions of the field and 
more precisely in the light-cone (LC) formalism because 
the Haag series can be used [16] and is well defined in 
terms of the product of <p(xi). In the full Euclidean 
metric there is no on-shell condition and <p(x) stays a 
D-dimensional Fourier transform with f(p 2 ) only. 
It might appear that there would be as many QFT's as 
eligible test functions. However some importants results 
from topological spaces analysis may be invoked. Some 
general definitions need first to be recalled : 

Definitions: 

a) An open covering of a topological space M is a 
family of open sets (?/»)»£/ such that M = Ujg/f/j. 

b) An open covering (?/»)»£/ is said to be locally finite 
if each point of M. has a neighbourhood which intersects 
only a finite number of the (Ui). 

c) A topological space M. is paracompact iff every 
open covering (t/j)j £ j admits a locally finite refinement 
(that is a locally finite open covering (Vj)j^j with the 
property that for all j £ J there exists i £ I such that 
Vj C 170- 

d) If (Ui)i£i is an open covering of M, a partition 
of unity subordinate to (t/»)»ei is a family ((3j)jej of 
positive continuous functions on M such that 

(i) For all X £ M. there exists a neighbourhood 
Ux of X such that all but a finite number of the (3j 
vanish on Ux- 

(ii) £ i€ i = ! for all X £ M. 

(iii) For all j £ J, there exists i £ I such that the 
closure of {X £ M : Pj(X) ^ 0} is contained in Ui. The 
closure of {X £ M : /3j(X) ^ 0} is called the support of 
ft. 

Remark: If and (flj)je.J are t w0 partitions of 

unity subordinate to (Ui)i^i and (Vj)j^j respectively, 
{ctii3j)(i.j)eixj is a partition of unity subordinate to 



(Ui n Vj)(ij)£i x j ■ Since this property is crucial in the 
sequel it is established in Appendix A. 

Then the following theorems hold: 

Theorems: 

A) Suppose M connected. Then M is paracompact 
iff the topology of A4 is metrisable iff (when M. is a 
diffcrentiable manifold) M. admits a Riemannian metric. 

B) If M. is a connected, paracompact differcntiablc 
manifold and (Ui) is an open covering, then there ex- 
ists a differentiable partition of unity subordinate to (Ui). 

Partitions of unity are useful because one can often use 
them to extend local constructions to the whole space. 
A well known case deals with the definition of integrals 
of differential forms [22]: let F be a differential form, 
((3i) is used to cut F into small pieces: Fi — [3iF and 
Y,F t = F. ¥ovVl£M one defines: J Q F := ^ J n . &F. 
(3i being of compact support each integral in the sum is 
finite and the overall result is independent of the choice 
of coordinates (atlas) on Oj and of the partition of unity. 

C) Localisation of distributions: J. Dicudonnc's GPT- 
theorem (Glucing-Picccs- Together) [18, 19] establishes 
the above properties for distribution functional. 

The metrisable Euclidean manifold being paracompact 
permits the use of partition of unity test functions in the 
definition of the field of Eq. (II. 5) . The resulting operator- 
valued functional is independent of the way this partition 
of unity is constructed. All constructions of the field 
with different partition of unity are therefore equivalent, 
thereby eliminating the initial arbitrariness in the choice 
of test functions. Then f(p 2 ) is 1 execpt in the bound- 
ary region where it is infinitely differcntiablc and goes to 
zero with all its derivatives. An immediate and impor- 
tant consequence of this construction is that all successive 
convolutions of the field ip(x) with the test function leave 
ip(x) unchanged. One has 



< T x ip,p>-- 



I 



d^-Vp f 2 (u 2 plP 2 ) 
(2tt)( d - 1 ) 2uj p 



a+e 



IpX 



+ a p e 



- vpx 1 



(II.6) 

Here / 2 stands for the product of two partitions of unity. 
According to the Remark above and Appendix A f 2 is 
also a partition of unity with the same support as / and 
ip(x) is not affected. 



B. Partition of unity and Poincare commutator 
algebra 

Without loss of generality, and to simplify the presen- 
tation, we consider here massive scalar fields at D = 4 
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and use the following conventions: dQ, k = m 
Vk 2 + m 2 ; [at, at] = {2i:f2uj k 8{k- k'). The energy- 



d 3 k 



;oj k = momentum tensor writes = d^ipd^p — ^^"{(dip) 2 — 



m 2 ip 2 ] and gives 



P° = H = J d£l k u k a+a % f 2 (w 2 ,k 2 ); P ] = \ j d£l k u> h h?[a+a % + a % a+]f 2 {oj 2 k , fc 2 ); (II.7) 
dn^a+ficolk^^-ia^lk 2 )); M jl =ij d£l k a+f{<4, k 2 )[k jw - fc,_](og/(^, fc 2 )). (H.8) 



r 



Evaluating the commutator of P J with y>(x) gives in turn 



(II.9) 

which is the expected result d^ip^x) because f 3 is also 
an equivalent partition of unity. For all other commuta- 
tors the reduction procedure due to the o's and a + 's is 
the usual one and since any power of the test function 
is an equivalent partition of unity the Poincare algebra 
is thereby satisfied. It is clear that any arbitrary test 
function will fail in that respect. 



III. BOGOLIUBOV-SHIRKOV-EPSTEIN- 
GLASER (BSEG) EXPANSION OF THE 
S-MATRIX AND LORENTZ INVARIANT 
EXTENSION OF SINGULAR DISTRIBUTIONS 

A. The BSEG procedure 

In the formal expansion of the S— matrix [2] the coeffi- 
cients are OPVD built out of products of free fields. Ac- 
cording to the above construction products of test func- 
tions appear then naturally and in the BSEG writing of 
S they are made explicit 

S (d) = 1 +^2^ fd i x 1 ..d i x n T n (x 1 , ..,x n )g{x 1 )..g{x n ). 

(mi) 

The test functions g(x) must have two properties: 

-i) vanish at infinite time in order to switch off the 
interaction. This is necessary to define asymptotically 
free states. It is a problem in equal time QFT because 
non-trivial vacua are excluded, but in LCQFT, because 
of the LC-vacuum structure, one can build a non 
perturbative theory on free field operators. 

-ii) restrict the propagation of fields to the interior of 
the light cone. 

Suppose that all time arguments of the set of vec- 
tors {x\,--- ,x m } are bigger than those of the set 
{x m+ i, • • • , x n } then the supports of the ensemble of test 



functions {g(x\), • • • , g(x m )} are all separated from those 
of {g(x m+ i, ■ ■ ■ ,g(x n )} . The requirement of causality 
leads to 

T n {x\ , • • • , x n ) — T m {x \ , , x m )T' n _ m (x m _)_i , • • • , x n ) . 

(III.2) 

As already mentioned in the introduction this condition 
completely fixes the dynamics of the system. If all ar- 
guments are different and ordered such that x\ > x\ > 
■ ■ ■ > by recursion it follows that T n (x\, ■ ■ ■ ,x n ) = 
T 1 (x 1 )T 1 {x 2 )---T 1 (x n ) and 



S(g)- 



00 1 r 

l+Y J -M A x 1 ..d A x n T[T 1 {x 1 )..T 1 (x n 

n=i nJ 



)]g{xi)..g{x n ). 



(III.3) 

It is important to note that the time ordering operation 
T cannot be made with ^-functions because the multi- 
plication of the two distributions 9(x) and Ti(x) at the 
same point is mathematically ill-defined. Epstein-Glaser 
analysis focusses then on the solution of this difficulty to 
which we now turn. 



Epstein-Glaser extension of singular 
distributions 



Because the functional integration of Eq.(II.l) may oc- 
cur in practice in a space of smaller dimension than the 
original dimension D we shall use d as a dimensional la- 
bel henceforth. Due to translation invariance the multi- 
plication problem just mentionned can be reduced to the 
study of distributions singular at the origin of ]R d . Let 
then f(X) be a C°° test function belonging to V(M d ) 
and T(X) a distribution belonging to V(M d \{0}) which 
we want to extend to the whole domain V(M d ). The 
singular order k of T(X) at the origin of 1R d is defined 
as 



k 



inf{s : lim X S T(XX) 



0}-d. 



(III.4) 



Epstein-Glaser original extension consists in performing 
an 'educated' Taylor surgery on the initial test function 
by throwing away the weigthed k-jet of f(X) at the ori- 
gin. Denoting by Rif(X) this Taylor's remainder, this 
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implies 

k 

R k f(X) = f(X)-f^ ]T ^-d a f(X)\ x=0 . (III.5) 

n=0 \ a\=n 

Here the notations are : a! = a 1 ! •••a!, \a\ — a 1 + 
■ ■ ■ + a d and d a = ■ ■ ■ df d . R k f(X) has the desired 
properties at the origin -it vanishes there to order k + 
1- but, due to uncontrolled UV-bchaviour, it does not 
belong to Schwartz-space V(1R d ). In order to be able 
to define a functional which can be used to perform a 
transposition in the sense of distributions Epstein and 
Glaser introduce a weight function w(X) belonging to 
V(R d ) with properties w(0) = 1, ^(O) = for 
<| a \< k. This allows to define a modified Taylor 
remainder eg: 

R k w f(x) = f(x)- W (x)j2 E ^r 9a /wi*= o - 

n=0|a|=n 

(III.6) 

Now the extension T W (X) of T(X) can be defined by 
transposition: 

<T w ,f>:=<T,R k w f>. (IIL7) 

The IR-extension T W (X) so obtained is not unique. We 
shall come back to this point . There is an abundant lit- 
erature on this original procedure of Epstein-Glaser and 
some important improvements were proposed recently 
[9, 10]. Our observation here [17] concerns the use of 
test functions belonging to T>(1R d ) and vanishing at the 
origin of 1R d with all their derivatives (dubbed super- 
regular or SRTF hereafter). In this case one has strictly 
f(X) = R k f(X), for all the derivative terms in Eq.(III.5) 
are zero. Then R k f(X) belongs also to T>(R d ), there is 
no need to introduce a weight function and the following 
set of equalities holds 

<TJ >=< T, R k f >=< R k T, f >=< f,f>. (III.8) 

According to Section II / stands here for a partition of 
unity which allows defining < T,f > from theorem C. 
Then, T becomes a regular distribution over the whole 
manifold on which / is just 1 everywhere. The way 
this is achieved in practice will be shown in Section 
IV. Important consequences follow from these relations. 
They have to do with the splitting procedure of dis- 
tributions into advanced and retarded parts, Lorentz 
covariance of the Epstein-Glaser analysis, the connection 
with BPHZ renormalization and consequently the analy- 
sis of the IR and UV behaviour of the undcrlcr lying QFT. 



C. Splitting into advanced and retarded parts of 
singular distributions 

As mentionned earlier the time ordering operation T 
in Eq.(III.3) cannot be made blindly without facing di- 



vergences related to ill-defined products of distributions. 
This ordering operation requires to split distributions T 
into advanced and retarded parts T a and T r respectively 
in the following way 

T{x\ , , x n ) — T r (x\ , • • • , x n ) T a {x \ , , x n ). 

(III.9) 

The support of T r (x\, • • • , x n ) (viz. T a (x\, • • • , x n ) )is 
the n-dimensional generalization of the closed foreward 
(backward) cone of {x\,--- ,x n }. Because of transla- 
tional invariance one can put x n = 0. If T(x\, • • • , x n ) 
were regular at Xi = 0,i = 1,- ■ ■ ,n~ 1 the splitting could 
be performed with 

n-l 

T r (x\ , .. , X n ) J_ \P ip^j *^n)-^('^l •En: ■• > ^n— 1 %n ; 0); 

i=l 
n-l 

T a (xi , .. , Xn) J^(*^n %j)T(xi x n7 .. , X n — i x n , 0) . 

i=l 

(111.10) 

On the contrary, if the product of ^-functions with the 
distribution T is ill-defined the splitting procedure has to 
be done with the help of Eq.(III.8). The defining equa- 
tions of T r and T a are then 

<T,6f> = < Of, I >=< T r , f >; (III.ll) 
<T,(l-0)f > = < {l-0)T,f>=<T a ,f>. 

T r and T a are simply obtained by multiplication of T with 
the corresponding ^-functions. With the original singular 
distribution T this would not have been possible. Finally 
the following identification is obtained 

f r = ef; T a = {l-e)T. (III. 12) 

It should be kept in mind that T r and T a are not unique, 
since, as stated above, T is not unique. 
If T(X) = T(x\,--- ,x n ) is causal -as is the case, if 
it is constructed with the BSEG procedure sketched 
in Section (III. A)- then the products of ^-functions 
in Eq.(III.lO) allow an essential simplification partic- 
ularly useful for calculations in momentum space [4]. 
There Scharf defines a 4(n — l)-dimcnsional vector 
v = (vi,V2,--- ,v n -\) consisting of n — 1 time- like 
four-vectors lying all in the interior of the foreward 
light-cone. The scalar product Vi ■ Xi = v®Xf — 1j-X 
is used to define the scalar v ■ X = YllZi v j ' Xj, where 
one puts X n = 0. Then v-X has the following properties: 

v ■ X > 0, \/Xj in the foreward light-cone F+ 

v ■ X < 0, \/Xj in the backward light-cone T~ 

Consequently v ■ X = defines an hyperplane that 
separates the causal support into advanced and retarded 
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parts. 

The following indcntities then hold: 

n-l 

6{v-x) = Y[e(x°-x° n ); XeT + 

i=i 

n-l 

6{vX) = Y[6(XZ-X$); XeT- (111.13) 
i=i 

These identities can only be used with distributions 
having causal supports. Otherwise the scalar products 
Vi ■ Xi would not have a unique sign related to advanced 
or retarded coordinates and v ■ X = could not be 
used to distinguish advanced and retarded parts of the 
supports. For later use we need the Fourier transform 
of 6(v ■ X), i.e. 6 v (p) where p = (p 1 ,p 2 ,--- ,Pn-i)- 
This multi-dimensional quantity can be reduced to a 
4-dimensional one by 

(1) choosing a coordinate system for which 
p = (pi,0, ■■■ ,0) where p\ — (p?,pi) (which can 
be obtained by an orthogonal transformation in 

(2) and choosing v = (v\, 0, • • • ,0); v\ = (l,~v) which 
yields 9(v ■ X) = 6(X° -xi-lf). Then 



UP) = (2^72 Jd D X ie ^e{Xl -%■-$) (111.14) 
Now with 

one obtains: 

T v (p) = (2^ D l^^—5^\^ -p°H). (111.15) 
p v + le 

The last equality means that ~p and ~~v come out paral- 
lel to each other. Of course physical results should not 
depend upon the choice of the arbitrary vector v. This 
will be verified in Section V. 



D. Lorentz invariance of Epstein-Glaser procedure 
and super-regular test functions. 

In the Taylor surgery above, under the action of an 
element A of the Lorentz group, derivatives transform as 

x a d a (Af) = /[A-^^/loA- 1 , 

= (A-^J^a^oA- 1 , (111.16) 

and in the general Taylor expansion x a d a (Af)(0) = 
(A~ 1 X) f3 dpf(0). Lorentz invariance is therefore violated 



in this procedure. However, as mentionned earlier, T(X) 
is determined up to a sum of derivatives of 6 functions 

^ (— 1) Q — j8^ a \X). But these 5-terms transform as 

\a\<k '" 

5( a >(AX) = [A- 1 *]^)^). Epstein-Glaser remedy 
is then to determine the a a 's to correct for the viola- 
tion due to derivatives. With SRTF vanishing at the 
origin with all their derivatives, because of the identity 
f(X) = Rif(X), Lorentz violating derivatives are just 
not there. Provided then that the space of test functions 
to be used is restricted to SRTF types, Lorentz invari- 
ance is satisfied from the outset. Formally T(X) is only 
determined up to the sum of derivatives of 5-terms but 
their contributions are redundant with super-regular test 
functions. Nevertheless its presence should not be forgot- 
ten and may prove instrumental in restoring some other 
broken symmetries [4] . 

E. BPHZ renormalization: a corollary of 
Epstein-Glaser procedure. 

One of the obstacles for a widespread application of 
the Epstein-Glaser approach was the non evident link 
to other renormalization schemes and in particular the 
question of its multipicative structure and its feasability 
in renormalization group studies. The chain of identities 
of Eq.(III.8) valid for SRTF allows to establish quite sim- 
ply the link of our approach with BPHZ renormalization. 
In the context of the Cesaro interpretation of diverging 
integrals [23] this link has been discussed by Prange [9] 
and Gracia-Bondia [10]. In the approach with SRTF the 
link to BPHZ is obtained as follows. 
Fourier and inverse Fourier transforms are defined as 

F(P) = Hf]( P ) = J J^^f(X); 

F(P) = T-\f]{p) = J 1 0^e i r x f(X). (111.17) 

Then, with /i a multi-index as defined after Eq.(III.5), 
the following relations hold 

0"/(O) = (-z)W(27r)- d / 2 < p^,¥ >, 
Jxi^ip) = (-t) M (2TT) d/2 5^(p). (111.18) 

The Fourier transform of the singular T(X) is well de- 
fined since the functional built with SRTF identical to 
their Taylor remainder is itself well defined. Therefore 
relations (III.8) and (III. 18) imply 

<tj> ee <T,i?f/>-<^(r),^- 1 (^/)> 

= <T{T) 1 R k 1 (T- l (f))> 

= <R k 1 f(T),f-\f) >, (111.19) 

that is 

T{f) = R\T(T). (111.20) 
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This is BPHZ substraction at zero momentum, up, as 
is well known, to an arbitrary polynomial in p originat- 
ing from the sum of (5-terms mentioned previously. It is 
known that for a non-massive QFT the BPHZ method 
faces infrared divergences. However, a mass scale can 
be introduced by doing substractions at some external 
momenta q ^ 0. For non-SRTF it amounts to choos- 
ing the weight function w(X) — e lqX . This is a legiti- 
mate choice provided the functional integral in Eq.(III.7) 
is given a meaning in the sense of Cesaro summability 
[10, 23]. However the situation is much simpler with 
SRTF. One has still e iqX f(X) = R k (e iqX f(X)) and the 
chain of equalities in Eq.(III.19) can be rewritten with 
this modification of /, leading now to BPHZ substrac- 
tion at momentum p = q. Thus BPHZ appears just as a 
special case of Epstein-Glaser method. In the next sec- 
tion we examine in more details UV and IR behaviours 
when using partition of unity test functions introduced 
in Section II. 



IV. IR AND UV EXTENSIONS OF T(X) WITH 
LAGRANGE'S FORMULA FOR TAYLOR'S 
REMAINDER AND PARTITION OF UNITY 

A. Lagrange's formula in the IR 

It is known from general functional analysis that Tay- 
lor's remainder can be expressed by Lagrange's formula. 

I 



Performing partial integrations in X, surface terms are 
avoided with SRTF. This is necessary to validate the 
usual operations with distributions. It turns out from 
the generic features of f(X) [17], detailed in Appendix 
B, that the t-integral has an effective lower cut-off at 
p,\\X\\ < 1 for jl < 1, where ||X|| < 1 is the norm of 
X. The extension is looked for in the IR region of the 
variable X. Since in the following X can be a spatial or 

I 



It is only recently [8] that its use was advocated in the 
Epstein-Glaser context of extension of singular distribu- 
tions and in the ensuing RG-analysis. With Lagrange's 
formula the transposition operation defining the exten- 
sion T(X) reduces to partial integrations in the func- 
tional integral. According to Lagrange R k f(X) can be 
written as 



vP r 1 

R k f(X) = (k + l) ]T hp / dt(l-t)% x) f(tX)]. 

(IV.l) 

For a SRTF <= V(M d ) it is easily verified that for any k 
Eq.(IV.l) expresses the identity f(X) = R k f(X). Since 
this identity is also valid for any power of f(X), the field 
construction and Poincare commutator algebra keep their 
initial properties with respect to the partition of unity 
(section II-(A,B)). One has then 



(IV.2) 

I 

a momentum variable "infrared" in the present context 
can mean UV (small distance) or IR (large distance) in 
physical terminology. Since for a momemtum variable 
the extension T < (X) is infrared and/or ultraviolet regu- 
lar by construction, the test function can now be replaced 
by 1 over the whole domain of integration [17]. One ob- 
tains: 



<TJ >=< T, R k J >= (k + 1) Yl J d d XT(X) [^ J dt(l t) k df tx) f{tX)] . 



\p\=k+i 



f<(X) = + 1) £ x [f f dt^T^)]. (IV.3) 

|/3|=fe+l ' J M\X\\ 



For an homogeneous distribution, that is T(j) = 
ti k + d )T(X), the i-integration can be carried out to give 
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~ X@ (—\ k 

T<(X) = (-) k (k + l) £ d^ x [—T(X)log^\\X\\)] + [ -^H k J2c^W^), 



(IV.4) 



r 



with H k 



C /3 



k 

E 

P =i 



(-i)'p +i > 



7 + + 1) an d 



J (||X||=1) T(X)X^. 
It is interesting to note how this result [24] parallels 
that of Ref.[10]. Here it is the behaviour of the SRTF 
at the origin which provides the lower bound in the 
t-intcgral and at the same time gives the identity 
f(X) = (1 - wiX^R^fiX) + w(X)R k f(X) Vw(X). 
It just corresponds to the T w operation of [10] applied 
on any test function <f>(X) G V(]R d ) with the choice 
w(X)=6(l-ii\\X\\). 

The power of the IR treatment through Eq.(IV.4) is 
shown in [10, 17] for the free massive scalar field the- 
ory developed in the (conventionally hopeless) mass per- 
turbation expansion. It leads through the exact resum- 
mation of the infinite mass-expansion series to the well 
known results given in terms of modified Bcssel functions 
./■^(toII V||). Another interesting test concerns integrable 
conformal field theories with additive interacting terms 
of massive dimension, which endure also conventionally 
untractable IR divergences if treated perturbatively [25]. 



Whatever the construction of the partition of unity 
underlying f(X) h is a parameter which may depend 
on \\X\\. The arguments underlying this observation 
are as follows. With h taken as a true constant the 
regularisation will be achieved in the usual cut-off like 
manner with well known symmetry preserving problems. 
In this case the final extension of the partition of unity 
to the whole manifold can only take place by finally 
letting h -the cut-off- going to infinity thereby losing 
the scaling information originally contained in the 
construction of the overlapping domains. Clearly using 
a partition of unity would be a useless artifact in this 
case. However with an ||V|| -dependant h, as we shall 
sec, the situation is very different: the final extension 
of the partition of unity to the whole manifold can be 
achieved altogether in a manner preserving symmetries 
and scaling informations. With an ||X|| -dependant h the 
consequences are: 



(j, 



-i)3 
2ll X\\ m 



| ^|| maa 

*9(\\x 



such that II X\ 



9{ 



\X\ 



i + MII*ll 



B. Extension of T(X) in the UV domain 

From the above IR analysis one could also explore the 
Fourier space UV domain. However the same formalism 
can be directly applied if T(X) gives rise to UV diver- 
gences in the absence of test functions. In QFT T(X) is 
in general a distribution in the variable |j X\\ only. Setting 
\\X\\ = X from now on, the domain of f(X) is the ball 
B 1+h (\\X\\) of radius 1 + h and f( a \l + h) = 0,Va > 0. 
With f(X) = f>(X), it also holds that 

f>{X) = -± ^ dt(l-t) k dl k+1 \t k f>(Xt)], 
I 



-ii) h > => M 2 ||^ll5(ll^ll) > 1 V ||X|| e 
[1.11*11 

max max )■> 



—Hi) from f > (tX) present in Lagrange's formula one 
ha8t<i±Mp=/xVlTO- 



In the definition of /i(||X||) a dimensionless scale factor 
/j? has been extracted from <7(||V||). Then one has 



<TJ>>= Jd d XT(X){-^^ ^(l-tM +1) [X k f>(Xt)]} 
= <T>,1>, 

I 



(IV.6) 



where in the last line the partial integrations in X have been performed. It gives the Z7U-extension T > (X) of 
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T(A) over the whole X domain on which / is just one 
everywhere. . 



T>(X): 



k\ 



X (k-d+l) d (k+l) [ X d T(X) 



^9(I1X||) 

dt 

T 



A direct application of this relation in the Euclidean met- 
ric is made for the scalar propagator A(x — y) which is 
diverging for D = 2 • • • 4 when x = y. These cases are im- 
portant to clarify the role of the presently unknown func- 
tion g(||A||), on the one hand in producing a regular UV- 



J 



extension T > (X) and, on the other hand, in the final RG- 
analysis with respect to the scale parameter /i present 
in /(A). Suppose now that T(X) ~ ||A||-( W+1 ) as 
||A|| -> oo, then / d^XT(X) diverges if (co + l-d) < 0. 
Thus k in Eq.(IVJ) should be such that k > (d- w - 1). 
Consider the Euclidean scalar propagator at x = y and 
for D = 4. The argument of the test function is dimen- 

2 

sionless. Hence we have X = j^, where A is an arbitrary 
scale which will prove irrelevant. Thus T(X) = XK } +m i 
and for D = 4 the space dimension in the A-variable is 
d = 2 and k — 1. We have then 



(p 2 + m 2 ) 



V,D=4 



X 2 



(XA 2 + m 2 ) 



dt 



r 



(IV.8) 



The full evaluation of Eq.(IV.8) is straightforeward. In 
keeping with the D = 2 case [17] the elimination of all 
non self-converging contributions in the A-intcgral re- 
quires that g(x) = a^" -1 ) (up to a multiplicative arbi- 
trary constant already taken into account as /j 2 ) and that 
the limit a — > 1_ is performed. Then h(x) = [i 2 x a — 1 
with < a < 1 is consistent with the generic construc- 
tion of f(X) (cf Appendix B). It implies also g(l) = 1 > 



J 



2^( 



(l-c) 



In the 



g(X max ) = 4r ie [i 2 > 1 and X max = (p 
limit a — > 1_ the upper integration limit in X extends 
then to infinity, the A— integral is self-converging and the 
test function can be taken to unity over the whole inte- 
gration domain. In this limit the propagator at x = y 
becomes 



A(0) = lim 



(47T) 



d A P f 2 ( P 2 ) 

(27r) 4 (p 2 +TO 2 ) 

> 2 - 1 - log( M 2 )), 



2m 4 



r 



d 4 p (p 2 - 1 - log( M 2 )) 
(2tt) 4 {p 2 + m 2 f 



(IV.9) 



which is the expected result with respect to the scale 
parameter [i [26]. One sees also that this parameter \i 
will always be present independantly of the infinitely 
many possible partitions of unity (cf section II) that can 
be used to build up /(A). In this sense fi is universal, 
hence its relevance for the final RG-analysis. The 
fundamental achievement of Lagrange's formula is that 
it perfoms the scaling analysis of the distribution T(A) 
through the integral on the t variable which incorporates 
the relevant scaling informations carried by /i. The 
analysis of cases involving loops (two- and four-points 
vertex functions) follows the same line and is detailed in 
Appendix C. 

Technically things are a bit more involved in 
Minkowski's metric. The test function becomes a 
function p(xo,~x) in coordinate space or /(p 2 ,,!? 2 ) in 
momentum space. As usual it is preferable to first 



carry out the integrations over xq or po respectively and 
then the remaining integrals. There are two possible 
situations: 

(1) The integrals over xq or pg are convergent. In this 
case the remaining integrations are as in the Euclidean 
space with dimension D — 1 and the extension is 
calculated correspondingly. As is shown in section V 
there is no divergence problem for the po-hrtegration in 
the calculation of the retarded and advanced extended 
distributions of causal distributions. 

(2) The integrals over xo or p diverge. This requires 
an extension of the distribution with respect to the 
dependance on xq or po . 

For practical reasons the calculations are usually done 
in momentum space and we restrict the discussion to this 
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case. In principle there can be an UV-divergence prob- 
lem of the po"i n tegration or there can be nonintegrable 
poles at finite values of po- The latter singularities are 
of IR-nature. Such situations can arise when extending 
non-causal distributions like powers of Feynman propa- 
gators. Once the po-extended distribution has been cal- 



culated the remaining integrations are to be done as in 

As a matter of illustration we recalculate the (Euclidean) 
result of eq.(IV.9) in Minkowski's metric (example of in- 
tegrable pole contribution). We want to give a meaning 
to the diverging integral A (0) at D = 2, 4 



J 



A(0) = 



d D p 



IHpIp 2 ) 



pi -p 2 



2ieu v 



r d( D -Vp r 

J 2u p J_ 



dp [ 

oo u p - po - ie 



+ ■ 



(IV.10) 



The po-integration cannot be done using contour inte- whole complex plane is not possible in general. However 
gration because the extension of the test function to the one can proceed as follows. One has 



PP 



dp 



Po ± 



dpo 



d Po} — x 1 



-Pa 



dpo 
= 0, 



dt 



f(p 2 t 2 , P 2 )} 



(IV.ll) 



r 



where a partial integration on p has been performed. 
Hence 



J — < 



dp 



P(pIp 2 ) 



±iirf 2 (u 2 ,p 2 ). (IV.12) 



po ± ui p T it 
At dimension D = 2 one obtains 

A( 0) = -nr r ^f(u 2 p ,p 2 ), 

J-oo Up 

= - in < — ./^wjj.p 2 ) >= -*7T < — 
Up ' Up 



(IV.13) 
1 > . 



Here d = 1, 



0, and fc = 0. Thus 



. 1 , d . dt 
(— ) = 



u v 



dp Ji t y/ p 2 + m 2 t 2 
1 1 



\Jp 2 + m 2 yjp 2 + m 2 ii 4 



(IV. 14) 



The end result for A(0) scales then as log[^ 2 ] as expected. 
The two contributions in Eq.(IV.14) may be recombincd 
to give a PV-type of substraction at the level of the prop- 



agator, as discussed more extensively in the next sub- 
section 



A(0) = / d^p^M^L, 
J p z — m A + le 



(IV.15) 



/ 



= d 2 p[ 



1 



p^ — ra A + le p A — m^/i 4 + it 
It is important to note that causality is restored since the 
non-causal diverging term in 8{x 2 ) now present in each 
contribution of Eq.(IV.15) just cancels out in the sub- 
straction. Here this is a consequence of "repairing" an 
ill defined Feynman propagator. In CPT it is a general 
feature that - by construction - diverging causality vio- 
lating terms are avoided right from the beginning [4] . At 
dimension D = 4 one has 



A(0) = -4i7r 2 A 3 / 
Jo 

= -4in 2 A 3 < 



XdXf 2 (X) 
y/X(XA 2 + rrPy' 



^X(XA 2 + m 2 ) 



1 > .(IV.16) 



Here d = 2, u — 0, and k = 1. Thus 



( 



y/X(XA 2 + m?) 



) = -d 2 x [X 



t 2 y/X(XA 2 + mH) 



(IV. 17) 
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The final t and X integrations give back the result of Setting \\X\ 
Eq.(IV.9). metric for D = 



= X as before, we have, in Euclidean 
2, d = 1 and k = which gives 



C. Transcription of Epstein-Glaser method in 
terms of Pauli-Villars type of substractions at the 
level of propagators 

Again we consider the scalar propagator A(x — y) in 
dimension D = 2 for the simplest demonstration of this 
transcription. It proceeds as follows. 
An alternative form of T y (X) is obtained through the 
change of variables Xt — > Y in the first line of Eq.(IV.6). 
It gives 



2 

f>(X) = tl)Hx^- d +^d^ +1) [x d [ dt 

k' Ji 



(1 



t (d+i) ^ t >\ 
(IV.18) 



(p 2 + m 2 ) 



alter 



= d x [xj^ 



dt 



1 



t (XA 2 + m 2 t) 1 
1 



(p 2 +m 2 ) (p 2 + m 2 n 2 ) 



This is a Pauli-Villars subtraction, however without 
any additionnal scalar field. The final momentum in- 
tegration gives the familiar RG-invariant result A(0) = 
j^pj \og(fi 2 ). The situation for D = 4 is a bit more intri- 
cate but also carries uscfull informations about the limi- 
tations of the P V-subtractions with respect to the scaling 
analysis. As seen above the dimension in the X variable 
is d = 2 and k = 1. The scalar propagator at x = y 
takes then a form somewhat different from that given in 
Eq.(IV.8), 



(p 2 + m 2 



alter 



t) 



1 



{XA 2 +m 2 t) 1 

1 2 ATA 2 



(IV.19) 



XA 2 + m 2 t (XA 2 + m 2 t) 2 
I 



+ 



X 2 A 2 



(XA 2 +m 2 t) 31 



It is verified that the final AT-integration gives the very 
same result for A(0) as in Eq.(IV.9). In Eq. (IV.19) the 
expression I(X, m) in brackets is regular at X = and 
behaves as ^ for X — > oo. Relation (IV.19) is a particu- 
lar type of PV substraction. Indeed it may be written as 
(ignoring t and A which are inessential for the argument) 



I{X,m) 



1 



■[1- 



yX + (i 1 X 2 + SX 



A 2 



X + m 2L X 

A + BX + CX 2 
(X + m 2 )(X + A 2 )(X + A 2 ) 



Al)(X + A 2 Y 
(IV.20) 



with a = 2,7 = l,/3 = 5 = e = and Ai = A2 = m. 
In the general form of Eq. (IV.20), imposing the fall off 
in for X — > 00 gives B = 0, C = 0. Without loss of 
generality one may take S = e = then (3 = A\ + (1 — 
a)A| , 7 = (a — 1) and A = A|(a — 1). Normalizing 



such that 



(m 2 -A^)(m 2 -A^) 

known result that a general rcgularisation of -5^ 

(A f 



1 determines a to give the 
may 



proceed through multiplication by II" =1 ( p2+A 2) 
decomposition in terms of PV contributions 



with the 



(A 2 - m 2 ) 



+ m 2 1=1 (p 2 +A 2 ) p 2 + 



Mj (A 2 -A 2 ) P 2 + A 2 ' 



(IV.21) 



r 



Under this form it is clear that each individual PV-term 
treated separctly would have to undergo Epstein-Glaser 
treatment to exhibit its scaling behaviour in terms of /1. 

It shows that even though F\(p 2 ) = ^2+™^ m ight 
be viewed as a special case of test function reducing to 
unity when {A^} — ► 00 and obeying Lagrange's formula 
Eq.(IV.l) with X = p 2 it does not have the intrinsic prop- 



erties of partition of unity test functions beeing essential 
for the scaling analysis as described in the preceeding 
paragraph. 
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V. CAUSAL SPLITTING OF DISTRIBUTIONS 
AND LINK WITH DISPERTION RELATION 
TECHNICS 

Here the calculations of the retarded extension of a 
singular distribution is worked out explicitely. The result 
shows an interesting link to substraction technics known 
from dispersion relations. 

The starting point is the following form of Lagrange's 
formula for SRTF's: 

0PO = (- + i)E i^r d^-^d^ixt)]. (v.i) 

|/3|=w+l f ' J0 



It yields T<(X) after partial integrations and taking into 
account the restriction brought about by (j){Xt) in the 
t— integral (cf Appendix B, Eq.(B.3) and the discussion 
thereafter) : 



t<(x) = (-!)(-+!) y d 



*fe^^)T(f)]. 



' t (d+w+l) ^ t 



(V.2) 



r 



It corresponds to Eq.(IV.3) but with a less refined lower 
bound on the t— integral. However, in keeping with the 
PV-type of substraction (c/Eqs. (IV. 14-15)), we shall see 
below that it provides also the interpretation of the sub- 



traction in the dispersion relations. The corresponding 
Fourier transform -which becomes UV— regulated [27], 
hence its notation- is: 



W(p) = ^±11 y \ p - [ 



f("+l) 



J dkO v (k)d%T(pt-k)]. 



(V.3) 



Using Eq.(III.15) for 9 v (k) this becomes: 



0). 



(V.4) 



With the change of variables p\t — = kf and after 
cj + 1 partial integrations in one obtains: 



2tt w! 



T(fcf,0, ••• ,0) 



(p? - fcf /i 2 + ie A i 2 )(-+ 1 )(p? - jfc* + ie) ' 



(V.5) 



where the final t— integration has been performed to give 
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hi. 



dt- 



(1-tf 



1 



lip? {p\t - kf + it)^+ 2 ) (oj + l)(p Q 1 -k' 1 + te)(p 1 -k' 1 fi 2 + ie^)^+ 1 )' 

I 



(V.6) 



Finally with the substitution k'i = spf in Eq.(V.5) and 
returning to the general frame [4] we get: 



— i (u 2 - 

T>(r>) = — — -1 



ds 



T{sp) 



(s - l//i 2 - it')^+ l ) (! - s + it') ' 
I 



(V.7) 



The retarded, extended distribution T>(p) satisfies 
an unsubstracted dispersion relation since T<(X) of 
Eq.(III.12) is a well-defined regular quantity. On the 
other hand the original retarded distribution T< = 9T 
is singular. A dispersion relation for this quantity would 
require w + 1 substractions. It is interesting to see that 
the explicit construction of the extended distribution in 



J 



terms of the original singular distribution leads to the fac- 
tor (s - 1/fi 2 - ie') (w+1) in the denominator of Eq.(V.7) 
which is characteristic of dispersion relations with u) + 1 
substractions, with one important difference: the sub- 
straction point is not arbitrary but as shown here- 
after. It is the scale /j, present in the SRTF which fixes 
this point. The calculation of the advanced, extended 

distribution T>(p) follows the same lines with a result 
similar to that of Eq.(V.7) with it — > —it. The difference 
of the retarded and advanced pieces reduces to 



-*=r -^=r U I 1 \ n ,,2 



n=0 



/oo 
dsT(sp)S (n) (s- 
-OO 



1 



t(p) -EE 



n=0|/3|= n 



\ g j TpW , P \ 

/3! V 2 



(V.; 



which is just the Taylor remainder of T(p — q) at 
9=4- -C? 5 (p) in the notation of [4])- that is the BPHZ 



substraction at q 



p 

7? 



discussed in section III-E. 



VI. NON-PERTURBATIVE LIGHT-FRONT QFT 
(LFQFT) AND THE BSEG PROCEDURE 

At various places in the text we commented on the 
necessity of using test functions when quantizing fields 
in Dirac's front form [28]. As stated in the introduction 
the restriction of quantum fields to a lightlike surface 
does not canonically exist and turns out to be the ma- 
jor mathematical problem to face first in dealing with 
LFQFT. The question was initially adressed in [13] and 
fully discussed in [14]. The outcome is the Fock expan- 
sion of the field operator written in LC- variables with test 



functions taking proper account of both the LC-induced 
IR singularity at k + = and of the UV behaviour in 
k + . For a scalar field theory the Haag series is well de- 
fined in terms of products of free fields and provides a 
non-perturbative LF-scheme to evaluate physical observ- 
ables by projection of the equations of motion and con- 
straints on different Fock sectors [16]. The procedure 
leads to coupled sets of integral equations for amplitudes 
in the Haag series, presenting non-standard renormaliza- 
tion difficulties common to all LF dynamics approaches 
[29, 30]. The main difficulty comes from the unavoid- 
able truncation of the coupled set of equations to a fi- 
nite number of amplitudes. This truncation jeopardizes 
the renormalization procedure unless special attention is 
given to the scaling behaviour of the neglected contribu- 
tions. It appears then that, if the BSEG procedure could 
be applied in the non-perturbative LFQFT context, then 
only finite quantities will occur in the course of calcula- 
tions and the amplitudes, instead of obeying coupled set 
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of equations, could be determined inductively. Thereby 
overlapping contributions and the renormalization diffi- 
culties mentioned above could be avoided. An immediate 
objection is that since the BSEG procedure is perturba- 
tive in essence it cannot fulfil our non-perturbative goals. 
However the same objection is also valid for the standard 



where H mt (x) is the interaction Hamiltonian and the 
time ordering is made with ^-functions, albeit the diver- 
gences occuring on the light cone when {xi-\ — Xi) 2 = 
[32] and subjects of Epstein-Glaser attention. The light 
front is defined by lu.x = 0, where u> is an arbitrary four 
vector such that lo 2 = 0. The expression of Eq.(VI.l) is 
then repesented in terms of the light front time a = uj.x 
for, if — Xi) 2 > 0, the signs of uj.{xi-\ — Xi) and 

— ti) are the same provided the time ordering is 
treated according to the generalized Epstein-Glaser con- 
siderations dcvclopped above which eliminates light-cone 
divergences. Hence H™ 1 = H mt , however with impor- 
tant caveats. We discuss first the case of interacting 
scalar fields. Obviously the asymptotic states have to be 
free bosonic fields ipo{x) as defined after Eq.(II.4). Since 
the LC-Lagrangian is singular, the quantization has to 
follow the specific rules for constrained systems [34, 35] . 
In the present case it means that in addition to the equa- 
tion of motion for the particle sector field (f(x) there ex- 
ists a constraint which lives only in the zero mode sector 
and is of nondynamical nature. If this constraint has a 
nontrivial solution ft then the total field <&(x) is a sum: 

= f(x) + 0. The zero mode field carries the infor- 
mation on nontrivial physics and is the LC-counterpart 
of a nontrivial vacuum in equal-time quantization. Con- 
sequently the field appearing in the interaction H mt (x) is 
the total field <I>(x). There is, however, a problem which 
cannot be overlooked: the zero mode Q following from 
the constraint can only be obtained approximately and 
recursively in an iterative procedure, which means that 
the interaction H (x) is a priori known only as far as f2 
is known. This was the bad news, shared with alterna- 
tive formulations such as that of Ref. [30] or the covariant 
light front dynamics approach of Ref. [29] ; the good news 
is that, if the coupling is large enough to support a zero 
mode - or in other words to generate a phase transition - 
already the first, algebraically very simple iteration step 
leads to a form for the zero mode giving a value for the 
critical coupling in good agreement with conventional, 
much more complex calculations and permitting to dis- 
cuss the critical behaviour near the phase transition [16]. 
Thus there is a good reason to hope that a small num- 
ber of iteration steps would lead to a proper description 



scheme with coupled equations and it is well known that 
the non-perturbative aspects of LFQFT are encoded in 
the constrained dynamics related to the singular nature 
of the light front Lagrangian. To be specific consider the 
standard expression for the S-matrix [31] 



(VI.1) 

I 

of non-perturbative physics. The presence of zero modes 
in the field operator does not spoil the causality issues of 
the theory: causality expressed through the Pauli- Jordan 
commutator function makes a statement about the mea- 
surability of fields at different space-time points. Since 
such a statement concerns only the particle sector of the 
theory, the zero mode contribution is not involved in the 
causal commutator. 

In fermionic theories the particle sector itself is split into 
a sum of independent and dependent degrees of freedom 
which leads -in gauge theories- to additional interactions. 
The zero mode problem concerns the gauge field itself and 
interaction terms built out of primary fermionic fields and 
behaving like scalar fields. 

The validation of these ideas in LFQFT gauge theories 
are undoubtly challenging. Already Yang-Mills theories 
have been investigated in some details within the stan- 
dard CPT framework [36] , paving the way to the imple- 
mentation of the CPT-LFQFT studies. They will be of 
considerable theoretical and even practical interests since 
the outcome is a rigourous non-perturbative alternative 
to lattice calculations, enjoying the numerous benefits of 
an LFQFT formulation put foreward for decades [37] . 



VII. CONCLUSIONS 

Our concerns in this paper were twofold. First to de- 
fine in a mathematically consistent way Quantum Fields 
as Operator Valued Distributions, focussing on the spe- 
cific test functions needed to achieve a generic QFT de- 
scription preserving basic Poincare and Lorentz invari- 
ances. In keeping with the general analysis of integrals 
of differential forms the necessary test functions take the 
form of partition of unity constructed from super regu- 
lar building blocks on compact subspaces. The second 
goal, related to the analysis of divergences crippling the 
usual pertubative QFT approach, was to show, following 
Epstein and Glaser analysis, the role of a proper treat- 
ment of causality in extensions of singular distributions. 
These extensions can be performed very elegantly with 
Lagrange's integral kernel representation combined with 
super regular test functions, for which the Lagrange re- 



S = l + Y.J Hr# mt 0d)^i - t 2 )H mt {x 2 ) ■ ■ ■ 9{t n - x - t n )H in \x n )d i x 1 ■ ■ ■ d 4 x n , 
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mainder of any order is identical to the test function it- 
self. The flexibility in the construction of the test func- 
tions allows to treat IR- and UV-divergent distributions 
with the same mathematical techniques. By the very na- 
ture of the construction an arbitrary scale, characterising 
the building blocks of the decomposition of unity, comes 
naturally into the picture. Such a scale is the corner stone 
permitting the use of extended distributions in renormal- 
ization group studies. Another important aspect of the 
decomposition of unity concept for the construction of 
the test functions is that it does not violate gauge invari- 
ance - contrary to other test function methods. A first 
treatment [38] in the QED context verifies this for the po- 
larisation tensor in dimension D = 2 and gives Fujikawa's 
analysis of the axial anomaly in dimensions D = 2,4 di- 
rectly from the presence of the test function as partiton 
of unity. A general analysis of gauge invariance is given 
in [4] and should be further developped in the present 
context. In coordinate space the extended distribution 



differs from the original singular distribution only on the 
support of the singularity whereas in momentum space 
it differs from the original one by subtractions at points 
determined by the intrinsic scale of the test function. 
Finally we presented arguments, based on the recent 
developments of Epstein-Glaser causal approach, which 
make it extremely plausible that a finite symmetry- 
preserving LCQFT could be envisaged on the basis of 
an iterative construction of the S'-matrix and a causality 
conditioned finite rcgularization using the OPVD treat- 
ments of fields. 
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Appendix 

APPENDIX A: COMPOSITION OF THE 
PRODUCT OF TWO PARTITIONS OF UNITY 

The purpose here is to give a proof of the Remark 
of section II, quoted quite frequently in the literature. 
Most of the time the proof is left to the reader but since 
the property is basic to section II it is useful to establish 
its validity. 

Property: If (ai)i e i and (Pj)jej are two partitions of 
unity subordinate to (Uk)keK and (Vi)i^l respectively, 
{otif3j)u t j)£ixJ is a partition of unity subordinate to 

1) Clearly (Uk H Vi)^j^ e KxL is a covering of M 

2) For X £ M one chooses a neighbourhood Ux (resp. 
Vx) of X such all but a finite number of on (resp. (3j) 
are zero on Ux (resp. Vx)- Then all but a finite number 
of the products oti(3j are zero on the neighbourhood 
Ux ("1 Vx de X. This prooves condition d(i) of section II 
for the family of functions (cKi/3j)(i,j)e/x J- 

3) Let ai, . . . , a. a et [3\, . . . ,[3b be the aj (resp. (3j 
which does not vanish on Ux (resp. Vx)- One has then 
for all x £ Ux, 

ai(x) + ... a A (x) = 1 
and for all x £ Vx, 

P 1 {x) + ...0 A (x) = 1 
which implies that for all x £ Ux H Vx, 

( ai (x) + ... a A {x))(l3i(x) + ... (3 A (x)) = 1. 



The products otipj appearing in this equation are the only 
one which does not vanish on C/^n Vx, which implies that 
for all x £ Ux n Vx, one has 

Y oti{x)Pj{x) = 1. 
(i,j)eixJ 

This proves condition d(ii)of section II for the family of 
functions {aif3j)(i,j) e i x j- 

4) Since the support of a product of fonction is the 
intersection of the supports, if for i £ I on chooses 
f(i) £ K such that suppcti C Ufu\ and for j £ J on 
chooses g(j) £ L such that suppfij C V g ^ it will hold 
that supp(ctil3j) C Uj{i) n V g (j). 

This proves condition d(iii)of section II for the family 
of functions (a»/3j)(i,j) e jx j- 

As a consequence the definition of the field given in 
Eq. (II-5) with / as a partition of unity is indcpendant of 
the way it is constructed since 

i^ 1 (*ij')€/xJ (j,i)£jxl je.J 



APPENDIX B: TEST FUNCTIONS AS 
PARTITION OF UNITY: CONSTRUCTION AND 
PROPERTIES. 

Since in the analysis of Section IV the generic test 
function depends only on one variable, the dimension- 
less variable \\X\\, here we only need to detail the sim- 
plest case of the decomposition of unity on the line. The 
family of functions {/3j} building up unity are based on 
an elementary function u with the property that for any 



17 



point x G [0, h] u(x) + u(h — x) = 1. This is the sim- 
plest realisation of the local finitcncss property of the 
{uj} family. By translation for any j G / and for any 
x G \jh,(j + l)h] onehasu(x — jh)+u((j + l)h — x)) = 1. 



Then (3j(x) = u(x - jh) and £jl-<x> 0j(x) = 1. The 
ovcrlaping subsets (Oj) covering the line are visualised 
in Fig.(Bl). 




Each individual /3j vanishes with all its derivatives out- 
side (Oj). Due to the independence of the OPVD func- 
tional on the construction of the partition of unity there 
is no need to consider more elaborate constructions in 
terms of the constituants functions {uj}. From the rela- 



jj(k) = dx u(x-jh) e~ lkx 

J(j-i)h 

2mr 



(B.l) 



= if k= — — , h if fc = 0, Vj. 
h 

This orthogonality property leads to the Fourier mode 
analysis in the distributional context. It is also interest- 
ing to note that the orthogonal Bloch functions approach 
used in the noncompact lattice formulation of gauge the- 
ories of Ref.[39] can be then rephrased in terms of the 
decompositon of unity, the lump functions there being 
particular realisations of the u(x — jh) with only a fi- 
nite number of derivatives vanishing outside (Oj). For 
completness we give explicitcly two constructions of the 
elementary function u. They are based on the regularis- 
ing function of L. Schwartz [18] 



p(x) = He 



if |a;| < 1, if \x\ > 1, (B.2) 



where M stems from normalising p(x) to unity. ui ,h \x) 
is then obtained from the convolution of the character- 



istic function {xj,h( x ) = 1 ^ jh < x < (j + l)h, 
elsewhere }, with the function p € (x) such that p € (x) = 
\p{-)- One has u^' h \x) ^= for x G [jh — e, (j + 
l)/i + £]andti^' i) (s) = l for x G [jh + e, (j + l)h - e]. 



tion u(x — jh) + u((j + l)h — xj) = 1 for x G [jh,{j + l)h] Th ^ n ^i 
it follows also that 



N-l (j,h) 



--0 



(x) is a partition of unity on = [0, L] 



with L — Nh. Another useful construction is achieved, 
for any v > 0, with 



u(x — jh) = < 



M I\x-jh\ dvex P 
for (j 



v u (h—v) u J ' 

l)h< x < (j + l)h; 



for \x — jh\ > h, 

(B.3) 

where M = J cfoexp[— v vh v y ]- It is verified that for 
any x G [jh, (j + l)h] one has indeed u{x — jh) + u((j + 
l)h — x)) = 1 and that (3j vanishes with all its derivatives 
outside (Oj) = [(j-l)h, (j + l)h}. 

In the scaled variable \\X\\ the final test function /(X) = 
f > (X) built in this way is then 



f(X) = f>(X) 



(l- X (||X|| + l,h)) for0< ||X|| <h 
1 for h< \\X\\ < 1 

X (||X||,h) forl<||X||<l- 
for || JSC || > 1 + 

(B.4) 



18 



f f 



1 




Fig.B2: The function f(X) from Eqs.(B4) and from u(x — jh) with v = 1 

a) left curves a — 0.95, dashed: /j 2 = 1.25; full line: ji 2 = 1.15; 

b) right curves fi 2 = 1.15, dashed: a = 0.95; full line: a = 0.85; 



Hence the function is just the very end part 

of the elementary function u building up unity, say 
w(||AT| — 1) in the last example. Since f(X) vanishes 
with all its derivatives at ||X|| = 1 + h, Lagrange's 
formula, Eq.(IV.5), is identically fulfilled. Moreover 
from its construction ,f{Xt) is different from zero 
if ||X||f < (1 + h). Hence the upper limit of the 
t-integral in Eq.(IV.5) is t max — ^pqp Due to the 
regularising exponential at the upper boundary it is 
easily seen on the two constructions presented above 
that f(X) still vanishes at an || V||-dependant upper 
boundary with all its derivatives and faster than any 
inverse power of \\X\\. This property is assumed to be 
preserved for the generic function x(||X||,/i(||X||)). One 
writes now h(X) = /j, 2 \\X\\g(\\X\\) — 1 where /i is an ar- 
bitrary dimensionless parameter extracted from g(||X||). 
It controls then the size of the overlapping regions for 
two consecutives (Oj).This is shown in Fig. (B2), where 
f(X), constructed from Eqs.(B4) and from u(x — jh) 
with h(X) = n 2 \\X\\ a — 1, is shown respectively at fixed 
a (viz. /it 2 ) for two values of fi 2 (viz.a). With the scale in 
the X variable used in Fig. (52) for < ||X|| < h f(X) 
rises so sharply to 1 that it is not distinguishable from a 
9— function. 

From the discussion in Section IV [i > 1 and 
tmax = = ^ 2 ||^!l ( " _1) > 1- The maximum 

value of ||X|| is (/i 2 )' 1 ^- 1 and tends to infinity when 
a tends to 1. For a distribution T(X) singular at 
the origin of IR d and with the representation of 
Eq.(B.3) the test function f < {X) behaves simply 
as w(X) = u(h - \\X\\) for \\X\\ < h. Then 
in the limit h — > w(X) tends to a step function 
9(\\X\\ - h). When a -> 1 h effectively tends to 
zero for values of II XII of the order of \ and then 



-(f) 



0( 



\x\\ 



9(t \\X\\(v 2 1)). 



Since the extension of T(X) is seeked in the IR and 
ft = ([i 2 — 1) > 0, jj,\\X\\ < 1 is a lower bound to the 
^-integral in Eq.(IV.3). 



APPENDIX C: EVALUATION OF 
LOOP-CONTRIBUTIONS WITH PARTITION OF 
UNITY TEST FUNCTIONS. 

We consider first the one-loop contribution I(k 2 ) to the 
four-point function of <E> 4 scalar field theory at D = 4. We 
shall work in Euclidean metric to simplify the derivation. 
I(k 2 ) reads 



f 2 (p 2 )f 2 ({p + k) 2 ) 



(C.l) 



)4 (p2 + m 2)((p + fc)2 + m 2y 

1 f d 4 q f 2 (( q + xk) 2 ))f 2 ((q-(l-x)k) 2 ) 
1 '"' '(27r) 4 (q 2 + k 2 x(l - x) + TO 2 ) 2 



The two cases k 2 = and k 2 ^ must be distinguished. 

p 2 

A? 

A 4 f°° XdXf 4 (X) 



In the first case, with X = one has 



7(0) 



(4tt) 2 J {Xh 2 + m 2 ) 2 



(C.2) 



With T(X) = {XA , X +m , r , VmT(X) - Y- Hence d 
1, oj = 0, k = and then 



T{X)=d x [- 



X 1 



2Xm 2 



L (XA 2 + m 2 ) 2S (VA 2 + m 2 ) 3 
It follows that 

2m 2 A 4 f°° XdX 



(C.3) 



7(0) = 



(4tt) 2 J (XA 2 + to 2 ) 3 J 1 t (Air) 



dt 



1 



log^ 2 ) 
. (C.4) 

In the second case, k 2 ^ 0, since the test function is unity 
for small values of its argument, in the UV regime x||fc|| 
and (1 — x) || fe|| can be disregarded with respect to ||g|| 
and the regulation of 7(fc 2 ) is provided by f 4 (q 2 ) only 
(This can be checked explicitely by Taylor expansion of 
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the f's, for in all explicitely finite contributions the limit 
a — * 1 can be taken, then / — > 1 everywhere and all 
terms involving derivatives of / with respect to g 2 do not 
contribute). Eq.(C.l) transforms to 



A 4 



(47r) 2 m 4 

A 4 
(47r) 2 m 4 



[dx [ 
Jo Jo 



dX- 



Xf\X) 



{X^ + ^x{l-x) + lf 
Y ™Y_f*[Y(£x(l-x) + l)] 



k 2 



dx 



(C.5) 



where, in the last equation, X has been changed to 

Y(^x(l -x) + 1). As before 

d = — 0, k = and I(k 2 ) reads now 



/(fc 2 ) = 



2A 4 



l dx I 
Jo Jo 



YdY 



+ I) 3 



(47r) 2 m 4 

/ -/ 4 [Ft( — za-aO + l)] (C.6) 
Ji t w 2 

As discussed in Appendix B the test function effectively 



cuts the t-integration such that t < —n — - 

6 - (-4-x(l-x) + l) 



the 



limit a — * 1, where Y max — > oo and the test function 
tends to unity over the whole integration domain in Y. 
Hence 



As expected this is another example that the infinite 
contributions in - present in the conventionnal formula- 
tion has been removed and replaced by the scale depen- 
dance originating from the overlapping domains building 
up the partition of unity. In a rcgularisation procedure 
by a cut-off the divergence is log(^y) and, according to 
Eq.(IV.15), the regulator mass is A = /j,m, hence the fi- 
nal result in log(yU 2 ). 

Instead of the derivation retained here one may use as 
well Schwinger's representation for the propagators. The 
test functions in this case just provide the necessary han- 
dling of divergences of the final integrals over Schwinger's 
parameter, with the very same result of Eq.(C.7).The 
way the regularisation works is most simply seen for A(0) 
at D = 4. One has 



A(0) 



(2ir) 4 q 2 
A 4 



(C.9) 



(47r) 2 m 2 



f°° , _„4 [°°du { 
Jo Vdye "Jo ^ 



Now for a SRTF it holds that 



= -i-{log( M 2 ) - f 1 dxlog[^x(l-x) + l}}. 
(47r) 2 J a m z 

^From the analysis in dimension D = 4 — e the result is 
known to be [7] 

I{k 2 ) = ^^-m- e/2 [dx[^x{l-x) + l]-^ (C.8) 
(4tt) 2 J m 2 



{An) 2 



(4 



Cte. 



f( y ~) 
u 



dt(l-t)d 2 t [tf( V -)] 

I u 

00 y(l 



(CIO) 



Hence the test function takes care of the divergence at 
u = 0, and after partial integrations all integrals are now 
finite giving the result of Eq. (IV. 9). 
We examplify now the general handling of loop contribu- 
tions with Schwinger's paramctrization for the two-point 
function T 2 [k] with two loops. It reads 



f Al A2 / 2 (g 2 )/ 2 (g 2 2 )/ 2 ((gl + g 2 - kf) 

J (2tt)8 (q 2 + m 2 ){ql + m 2 ){{ qi +q2-k) 2 + m 2 ) 

f Al A2 A 3 6(q 3 + k - q x - g 2 )/ 2 (g 2 )/ 2 (g 2 2 )/ 2 (g 2 ) 
J (2tt) 8 (g 2 +m 2 )(g 2 2 + m 2 )(g 2 +m 2 ) 

-J dye l\J (2^ (qf + m 2 ) 

= J d 4 ye* k -yf[ J da ie - a ^ 2 J ^ e -<*«?-*«-» f(qf). (C.ll) 



In the last relation / 2 (g 2 ) depends in fact only upon the argument X { = jfe. Then with the change of integration 
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variable onq% — > <^ 2 the argument of the test function 
becomes ^ A2 which implies (c/ Appendix B) ctjA 2 > 
^j-. The integrals over ctj have a lower boundary cut- 
off at ^p-, just as in the conventionnal treatment [26]. 
However using instead Lagrange's formula for the test 
functions gives a finite result, whose dependence on the 
scale parameter \x can be inferred from the conventionnal 
analysis with the replacement A — > p,m. 



APPENDIX D: RECOLLECTION OF 
LAGRANGE'S TYPES OF INTEGRAL 
FORMULAE FOR SUPER-REGULAR TEST 
FUNCTIONS. 



In this section we collect integral formulae used alter- 
natively in the UV and IR regimes at different places in 
the main text and show their equivalence. 
The first one encountered in the IR is given in Eq.(IV.l). 
At D — 1, after k partial integrations, the following iden- 
tity also holds for SRTF 

1 f°° d k+1 

poo 

= -- jr I dt(l 



Transcribed to arbitrary dimension and with X = p this 
takes the form 




(D.l) 



It is easy to check that an equivalent representation of 
f(p) is 



/>(p )s - (fc+ i) ^ r*!^^^/^)]- 

\f3\=k+l J1 ' H ' 

(D.2) 

Here d = D — 1. The same restriction as discussed ear- 
lier applies to the t-integral. After partial integration in 
the UV-extension of T(X) and in the limit a — > 1 the 
upper limit in the t- integration is just [i 2 . The Fourier 
transform <f>(x) to the conventionnal x-space representa- 
tion gives an alternative interpretation of Lagrange's for- 
mula, Eq.(IV.l), valid in the IR. Normalising the Fourier- 
integral with (27r)( D / 2 ) and after the change of integra- 
tion variable t — > 1/t it reads 



|/9|=fc+l P ' J ° 

(D.3) 

which is Eq.(IV.l). However when used in the analysis of 
a distribution singular in the IR and in the limit a — > 1 
the lower limit in the i-integration comes from the upper 
limit in the original p- integral and is now just With 
this interpretation Eq.(IV.l) is particularly useful in es- 
tablishing generalized dispersion-type relations for causal 
distributions in relation with the findings of [4] . 



